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I. INTRODUCTION 



Bose-Einstein condensates (BECs) are a promising 
physical system with which to realize the analogue grav- 
ity program L l|. This field advocates the observation of 
the analogues of phenomena such as: Hawking radia- 
tion cosmological particle production 0] and super- 
radiance [1]. The analogies follow from the equivalence 
between the equations of motion for a scalar quantum 
field in curved space-time and for the phonons of a flow- 
ing BEC in the hydrodynamic regime [2]. The interest 
in the analogue program arises because of the lack of 
prospects for experiments in the gravitational domain. 

A sonic horizon in a BEC is a surface on which the nor- 
mal component of the condensate flow exceeds the speed 
of sound. It is analogous to the event horizon of a black 
hole. Sonic horizons are predicted to Hawking radiate 
phonons Q. A variety of methods to create sonic hori- 
zons in BECs have been suggested, most of which gener- 
ate both black and white hole horizons |1, [f| 0, H, 111 • A 
white hole horizon is a surface on which the condensate 
flow drops back below the speed of sound. Due to the in- 
stabilities associated with white holes, methods which do 
not generate them are attractive, such as outflow from a 
reservoir [lfj| . In this paper we study a simplified version 
of the latter method, under realistic experimental condi- 
tions, including three-dimensions and three-body loss. 

In general relativity, Hawking radiation [lj, QJI is a 
quantum effect whereby a black hole horizon emits a 
thermal spectrum of particles. In the analogous sonic 
horizon case the emission temperature is given, for one- 
dimensional flow, by Q 
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v is the flow speed, c the speed of sound, M = v/c the 
Mach number and Xh denotes the position of the horizon, 
where v — c. 

We find that three-body loss severely limits the Hawk- 
ing temperature that can be achieved in BECs, indepen- 
dent of the detailed implementation of the sonic hori- 
zon. Amongst helium, rubidium, cesium and sodium, 
the latter allows the highest Hawking temperature, for 
fixed total atom number loss, because of its relatively 
low three-body loss rate and high two-body interaction 
strength. 

This paper is organized as follows. In section [Til we 
introduce our notation and formalism, after which we de- 
rive a general limit on the analogue Hawking temperature 
in section Hill Numerical studies of horizons generated by 
reservoir outflow are presented in section Hvl 



II. METHODOLOGY 

The dynamics of the BEC wavefunction cf>, including 
three-body losses, is described by the Gross-Pitaevskii 
equation (GPE) 
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Here m is the atomic mass, U the interaction related to 
the scattering length a s by U = 4wh 2 a s /m and K3 the 
three-body loss coefficient. V denotes an external poten- 
tial. Using 4> — yfn exp (i-d) , we can extract the atom 
number density n and velocity v = hS/d/m. Eq. ([3|) 
can be rewritten in terms of n and v, resulting in con- 
tinuity and Bernoulli equations for the BEC. For now 
we assume K 3 — and that the condensate does not 
vary on length scales shorter than the healing length 
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£ = l/^f%ima s = K/{y2mc). We also assume that the 
system can be modelled in one dimension: for example, 
flow in a narrow confining potential with longitudinal co- 
ordinate x and effective cross-section A{x). We can then 
write the approximate equations (To| : 



J = n(x)v(x)A(x), 

/! = -mv 2 (i) + Un(x) + V(x). 



(4) 
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Here /i is the chemical potential, J the particle flux and 
v the flow speed v = |v|. The local speed of sound c = 
y/Un/m. Without assuming the absence of short length 
scale variations in the condensate variables, the so called 
quantum pressure Q(x) — V y^" would contribute 
to the rhs. of Eq. ((5|). Combining Eqs. and (JS}, we 
obtain [Icl |: 
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The function q(x) contains all information about the 
trapping configuration. The quantity s(x) equals the 
speed of sound at the horizon s(xh) = c(xh), see Eq. Q. 
Differentiating with respect to x, solving for the partial 
derivative with respect to position M'(x), and applying 
L'Hospital's rule, we obtain for the horizon temperature, 
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where all variables are evaluated at the horizon. To reach 
this simple form we used M(xh) — 1, q(xh) — 3/2, see 
Eq. ([6]), and q'(xh} = [10]. The latter implies the reg- 
ularity condition [lj] 



{e 2 A'/A - V'/m) \ x = Xh = 0. 
From Eq. (J7J) we get 



!=■ (log A) 



4 ,V 
--(logA)'— 



(9) 



(10) 



While analytic expressions for the horizon temperature 
in the cases of constant potential V, or constant area A, 
have been previously obtained 0, [l3] , Eqs. §8§ and (fTU|) 
are valid where both vary. 

In several simulations we use the one dimensional 
Gross-Pitaevskii equation. For these cases, we assume 
that the BEC is harmonically confined in the transverse 
directions, with a large trapping frequency lu±. We 
can then perform a dimensional reduction of Eq. ((3j, 
assuming a fixed Gaussian transverse wave function, 
<j>(x,r,t) = exp(-r 2 /(2a 2 osc ))/ x /7:a 2 osc ip(x), where r is 
the transverse radial coordinate and a osc = \J h/(mu}±_). 
Then, we obtain an effective one-dimensional interac- 
tion strength Um = U/(2-Ka 2 osc ) and loss coefficient 
K 3AD =K 3 /(3n 2 a 4 osc ). 



All numerical solutions were obtained using dimension- 
less units (tilded in the following) , with </> = <f> a% 2 , t — 
tu±, x = x/a osc . Then U = 4%a s /a osc , U 1D = U/(2%), 
K 3 = K 3 /(LJ ± a 6 osc ) and K 3 , 1D = K 3 /(3n 2 ). The explicit 
form of the one-dimensional GPE after conversion to di- 
mensionless units is 
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Note that the ID density n\o - 
n = 1 0| 2 are connected by hiD 
dimcnsionless Hawking temperature 
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satisfying T H = u> ± hf H /k B = oj±7.6 x 10~ 12 f H . When 
presenting our results using dimensionless units, we al- 
ways state a concrete atomic system to which they per- 
tain. The purpose of the dimensionless units is then to 
allow simple conversion of the results between different 
atomic species, which is possible for our results without 
three-body losses. 

As our focus is on experimental realism, we base our 
results on numerical solutions of the full GPE. However, 
we make use of the approximations Eqs. (|4])- (ll0p to guide 
these numerical simulations. In particular, we can use 
Eq. ijHJ) to estimate the analogue Hawking temperature 
from a given potential, and the regularity condition 
to predict the location of the sonic horizon. 



III. LIMITS ON THE ANALOGUE HAWKING 
TEMPERATURE 

The analogy between quantum fields in curved space- 
time and the excitations in a BEC requires both the 
condensate and excitations to be in the hydrodynamic 
15]. The former requires variations in the den- 



regime 

sity and speed to occur only on length scales much larger 
than the healing length £, and the latter confines our 
analysis to phononic modes, i.e. those whose wave num- 
bers q fulfill q£ < 1. 

Here, we show that for a condensate horizon in the 
hydrodynamic regime, the achievable analogue Hawking 
temperature is limited by three-body loss processes. We 
formalize the condition that the flow has to remain in the 
hydrodynamic regime by demanding 
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for D ^> 1, where D is a constant chosen to specify the 
flow speed gradient we are judging to be acceptable. If 
the Mach number profile is roughly linear near the hori- 
zon, we are requiring that M varies at most by the frac- 
tion 1/D within the space of one healing length £. In- 
serting Eq. (fT3"|) and the expression for £ into Eq. (fTJ) we 
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obtain 



Un 



(14) 



All quantities are evaluated at the horizon. Hence high 
analogue Hawking temperatures require dense conden- 
sates. However, losses limit the achievable densities. 

To determine a quantitative limit on the analogue 
Hawking temperature, we investigate the effect of the 
loss in detail. We consider a flow that is harmonically 
confined with frequency uj± in the transverse directions, 
and is in either of two regimes: We call the condensate 
quasi-one dimensional (Q1D) when fi <C frw±. In this 
case transverse excitations are "frozen" , and we can fac- 
tor the BEC wave function into longitudinal and trans- 
verse parts; the latter remaining in the oscillator ground 
state [16| . The velocity is then exclusively in the longitu- 
dinal x direction. Our second regime is fi which 
we call "transverse Thomas- Fermi" (TTF) . Then we can 
consider the flow in the Thomas-Fermi approximation 
[ItJ and the velocity can in general have a transverse 
component. We assume the density profile 



(15) 



n q id(x, r, t) = n (x, t) exp - 



in the Q1D case. Note that, in this section only, the 
normalisation of the transverse wave function differs from 
that presented in section [II] to give no(x,t) the meaning 
of the actual time dependent peak density. 
In the TTF regime we assume 
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with n tt { (x, r, t) = 0, where Eq. (TIT))) would be negative. 

If we assume that the transverse structure of the con- 
densate is always well described by either of the profiles 
defined above, we can integrate out the transverse coor- 
dinates. In this case, the effective loss-equation for the 
peak density is 
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with «3 = 3 in the Q1D and a 3 = 4 in the TTF regime. 

For short times or weak loss, we can approximate 
Eq. (Ti"T|) by a finite difference equation. We now de- 
mand that losses at most reduce the peak density by a 
fraction / of the initial value within a time interval At. 
This corresponds to 



^n 3 At < /n . 
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We hence define the maximal allowed peak density under 
these conditions: 
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Finally we substitute Eq. (fT9j) for the density into 
Eq. fll]) and obtain 
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This limit on the analogue Hawking temperature is one 
of our key results. 



A. Comparison of BEC- Atom Species 

In the light of the preceding result, we now compare 
atomic species commonly employed in BEC experiments: 
4 He, 23 Na, 87 Rb and 137 Cs. The relevant parameters for 
these atoms are listed in table [I] 



atom 


U x 10 50 [Jm 3 ] 


K 2 x 10 20 [m 3 /s] 


Ks x 10 42 [m 6 /s] 


4 He 


15.7 


2, Ref. [18] 


9000, Ref. [18] 


23 Na 


1 


0.053, Ref. [19] 


2.12, Ref. [19] 


87 Rb 


0.5 


150, Ref. [20J 


32, Ref. [21] 


137 Cs 


0.66 


small 


130, Ref. [22] 



TABLE I: Parameters for atomic species. K2 is the two-body 
loss coeffcient. The values for cesium assume a magnetic field 
of 23 G, where the ratio of scattering length to three-body 
loss is most favorable 12211. 



According to Eq. ([20]) the combination of atomic pa- 
rameters that determines the limiting Hawking tempera- 
ture is U/^/Kl. For 4 Hc, 23 Na, 87 Rb and 137 Cs this has 
the respective values: 1.7, 6.9, 0.9, and 0.6, in units of 
10~ 30 Js 1/2 . The limiting Hawking temperature is four 
to ten times higher for 23 Na than for the other species. 

The scattering length a s , and hence the interaction 
strength U, can be tuned towards higher, more favorable 
values with the use of Feshbach resonances. However in 
the large a s regime, the three-body loss rate increases as 
a 4 . [22J . According to Eq. (f2"0"|) this rules out an improve- 
ment of Th using Feshbach resonances. 

In order to obtain numerical temperatures from 
Eq. (|20p we must specify the maximum Mach number 
gradient and the maximum acceptable loss rate con- 
straints, which are determined by the parameters D and 
f / At respectively. We choose D = 20, which places 
the sonic horizon marginally within the hydrodynamic 
regime. Note that to use the eikonal approximation for a 
phonon excitation near the horizon [23j, the condensate 
flow should vary little within one wavelength A = 2ir/q, 
where q is the wave number. For the shortest wavelength 
phonons l/q ~ £. Hence we require small Mach num- 
ber variations within A ~ 2-7r£ which is barely true for 
D = 20 ~ 67r. The necessity to consider phononic modes 
places stronger constraints on the flow variation than the 
requirement that the flow itself be in the hydrodynamic 
regime. Eq. (|20j) is weakly dependent on the maximum 
loss rate f /At. We choose to consider a loss of 10% of 
the peak density (/ = 0.1) within a time At = 50 ms, 
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so that //At = 2 s 1 . The maximum Hawking temper- 
atures using these criteria are listed in table [UJ 



atom 


n max x 10" 19 [m" 3 ] 


£htm] 


T„ [n 


K] T„/T crit 


4 He 


3.0 


0.42 


3.8 


0.10% 


23 Na 


194 


0.086 


16 


0.15% 


87 Rb 


50 


0.12 


2.2 


0.19% 


137 Cs 


25 


0.12 


1.4 


0.30% 



TABLE II: Achievable Hawking temperatures according to 
Eq. pQ|l. calculated using q 3 = 4, At = 0.05 s, / = 0.1 and 
D — 20. Also shown are the condensate densities and as- 
sociated healing lengths allowed by criterion Eq. (|19[l . The 
temperatures are also given in relation to the critical temper- 
ature T crit = (n/C(3/2)) 3/2 2irh 2 /(mk B ). 

The values vary by more than an order of magnitude 
between the atomic species, with sodium exhibiting the 
highest temperature of 16 nK. We also give Tjj in relation 
to the condensation temperature [HJ|. For the densities 
in column one of table [TTJ two-body losses are negligible 
compared to three-body losses for the values of K2 in 
table UJ 

The temperatures in table |TTJ compare well with the 
15 nK given for sodium in Ref. [lOj]. The values are low 
in spite of the optimistic values for At, D and / that 
we have used. By allowing large losses (/ = 0.1) in a 
short time (At = 50 ms), we are assuming quick phonon 
detection. Slower detection, and moving the horizon 
safely into the hydrodynamic regime, would require all 
three values to be adjusted in the direction producing de- 
creased Hawking temperature. Hence fast phonon detec- 
tion schemes such as proposed in Ref. [25| are attractive. 
Another useful tool would be suppression of inelastic col- 
lisions [2^ | . 

An interesting direction for future research is moti- 
vated by the low temperatures of table [TTJ Determining 
whether an emission of quanta also takes place in con- 
densates in which exceeds the hydrodynamic bound 
of Eq. (fT"3| (or D ~ 1) and how visible this radiation 
would be. A candidate for such a condensate flow is that 
past a gray soliton Q. As for D ~ 1 strong deviations 
from a thermal spectrum are expected outside of the hy- 
drodynamic regime, the use of Eq. JIJ for the radiation 
temperature no longer makes sense there. 

It has also been suggested to realise the analogue 
Hawking effect in degenerate Fermi gases [27], HH, for 
which three-body recombination would be suppressed. 

IV. RESERVOIR OUTFLOW 

Many horizon systems suggested in the literature pos- 
sess a white hole in addition to the black hole @, [E 0, H, Hi 
. A white hole is the surface where the flow speed returns 
below the speed of sound. It has been shown that the con- 
densate develops dynamical instabilities in the presence 
of a white hole horizon 0, [H, HO] ■ While these can be 



suppressed by appropriate boundary conditions (30| . in 
a realistic experiment they are likely to be present. In a 
range of numerical investigations of sonic horizons in ID 
and 2D, we found dynamical instability in all cases con- 
taining a white hole. It manifests itself in the repeated 
creation of gray solitons [9( . 

The reservoir outflow scenario of Giovanazzi et al. [To( | 
has a black hole horizon without a white hole horizon. 
They reported, using one-dimensional simulations with- 
out losses, a setup with dynamically stable horizon for- 
mation. We analyze a simplified version of this model 
in greater detail, as it is the most promising method for 
experimentally realizing a stable sonic horizon. We first 
outline the physical setup and then investigate details 
of the proposed experiment. Finally, we will put all the 
pieces together and present a complete 3D simulation 
covering all aspects of the outflow scenario. 



A. Outflow Scenario 

Consider a BEC reservoir from which the condensate 
is allowed to flow out. For example, consider the atoms 
confined between two potential humps in one dimension. 
If one of these humps is sufficiently weak, the condensate 
can leak out, with the BEC becoming supersonic and re- 
maining so, avoiding a white hole. We develop this idea 
of Giovanazzi et al. [10( in several ways. First, we simu- 
late the horizon in 3D, exploiting cylindrical symmetry. 
In particular we do not always assume a quasi- ID sit- 
uation, and hence study the transverse structure of the 
horizon. Second, as our discussion in section IIIII showed 
the importance of three-body losses, we include them in 
our simulations. Finally, we evaluate whether the optical 
piston described in Ref. [Io| is required in practice. 

A sketch of the outflow scenario is given in Fig. [T] 
The condensate is initially confined in an elongated, cigar 
shaped trap, sliced by blue detuned laser sheets which act 
as endcaps. This defines our reservoir. The sheets super- 
impose to achieve approximately top-hat shaped poten- 
tials. The height of this superposition potential V^ up is 
larger than the bulk condensate interaction energy Un^g , 
hence the atoms are confined, see Eq. ([5]). n^g is the 
background density within the reservoir, i.e. where the 
potential vanishes. To initiate outflow, all but one of the 
right-hand sheets are suddenly removed. The height of 
the remaining single sheet Vq is less than that required 
for confinement so that the BEC starts to stream out. 
The purpose of the superposition method is to avoid a 
rapid temporal change in the potential within the region 
where the condensate density is significant, which would 
disrupt the reservoir. 

In an experiment one could just let the condensate 
propagate supersonically along the waveguide created 
by the transverse confinement. However, the numerical 
treatment requires absorption at the grid edge, for which 
we employ a simple imaginary potential. 

Taking all these elements together, the potential that 
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FIG. 1: (Color online) (a) Sketch of an elongated condensate 
(gray), sliced by multiple blue detuned laser sheets (blue), (b) 
These sheets (blue) add up to form barriers (black, dotted), 
which are impenetrable for the condensate. The red, dashed 
line is the interaction energy Uh. The left sheets act as a 
permanent endcap. To initiate supersonic outflow from this 
reservoir, the optical potentials sketched in light blue in panels 
(a) and (b) are suddenly removed. 



we apply for our dynamics is 
1 Nl 

V(f,x,t) = -r 2 + ^2v G (y ,xi(t) - nd,a) 



a hump potential just high enough to confine the BEC 
by itself, Vq = ft. Then, simultaneously with the change 
in potential, we use the left sheets as a moving "optical 
piston" : We let xi — xi(t) — xi(0) + v p t. The sheets move 
with the velocity v p to the right and press the condensate 
out of the reservoir. When we model this, we additio nally 
phase imprint the velocity v p onto the condensate [3l[ , 
which makes the outflow smoother. 

When we consider three-dimensional situations, we 
solve the dimensionless variant of Eq. §3§ using cylin- 
drical coordinates f, x, ip. Besides small details in the 
shape of the laser sheets, which can be seen in Fig.[T](a), 
our problem is cylindrically symmetric. We neglect these 
details and assume independence of the potentials and 
wave function of if. The simulations are thus done on 
a two dimensional grid. Nonetheless we refer to them 
as "3D simulations", as their results should pertain to 
the fully three-dimensional scenario due to the approx- 
imate symmetry. We find the condensate initial state 
using imaginary time evolution 32]. The condensate is 
assumed to be in the ground state before t = 0. 

If we consider the effective cross-sectional area A(x) to 
be constant, the dimensionless analogue Hawking tem- 
perature at the hump potential (for t > 0) can be readily 
evaluated from Eqs. (T]), © and (EH). One finds [l(J: 



Z7T 



3 d 2 V 
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We know from section Mil that Hawking temperatures 
increase with density. This seems not to be the case in 
Eq. (|22[) . however there is a disguised dependence on n 
as we are considering the case where Vq j» fi = Un b g. 



B. Speed of Sound at the Horizon 



+ V G {V ,x ri a) + Vg(- H,x oc ,a- c) 



N r 



6{-t)J2VG(V ,x r +nd,a), 



n=l 

V G (W,xo,* ) = Wexp(-2(£-x ) 2 /a 2 ). (21) 

Vo denotes the amplitude and a the width of an individ- 
ual Gaussian sheet, xiit) is the position of the innermost 
sheet in the left barrier and x r the same on the right, d is 
the separation between the sheets. N t / r are the numbers 
of sheets in the left and right barriers. Usually xi(t) = xi 
is constant. The absorption is parametrized by 7, x oc 
and a oc . The potential Vg(Vq, x r , a) is distinct from the 
other members of the right top-hat array, as it is the 
one that remains after t = 0. We will refer to this sheet 
as the "hump potential". 9(— t) denotes the Heaviside 
step function. For our one dimensional simulations, the 
potential is given by Eq. (j2T|) with f = 0. 

Instead of initiating outflow by lowering the barrier on 
the right, we can also proceed as in Ref. [10(. We then use 



To obtain the most general temperature limit in sec- 
tion IIIIl we have previously derived the speed of sound 
at the horizon directly from the maximal density allowed 
by three body losses. In the outflow scenario presented 
in this section this overestimates the limit, as the density 
is usually much lower at the horizon than in the bulk 
condensate. Specifically we can evaluate Eq. © at the 
horizon: 



3 3 
H = -mc(x h ) 2 + V(x h ) = -mc{x h ) 2 + V 



(23) 



In the first step we have used v = c = yJUn/m at the 
horizon and in the second step the fact that the horizon 
is always established at the peak of the hump potential 
[l0|] . By the same equation ([5]) we can estimate the chem- 
ical potential as /i = XJn\, g since the velocity of the the 
bulk condensate is small and the potential vanishes there. 
Solving Eq. ([23]) for c(xh) 2 we obtain from Eq. ([14)): 



V2(Un bg -V ) 
37rfcs-D 



(24) 
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Hence the strength of the hump potential Vq has to be 
carefully adjusted. It should be large to increase the 
Hawking temperature Eq. (|22| but must be substantially 
less than the bulk pressure U nb g for the temperature not 
to violate Eq. (|24|) . Our numerical results in the follow- 
ing sections are chosen near these limits: We consider 
bulk densities as large as allowed by losses with a strong 
hump potential. As a consequence the flow at the hori- 
zon will show small deviations from the hydrodynamic 
regime. We provide the healing length in the background 
condensate the healing length at the horizon £/j as 
well as the maximal temperature T* allowed by Eq. (TT4|) 
(for D = 1) to relate the horizons to the hydrodynamic 



conditions. In dimensionless units: T* 



~;2 



(x h )/{V2n). 



Outflow Initiation 



In this section we compare the two different methods 
for initiating the outflow from the reservoir: (i) A sud- 
den change of the right barrier potential from a maxi- 
mal strength of V sup > fiiD to Vq < jliD- (h) A sud- 
den change of the right barrier potential from a maximal 
strength of V^ up > fim to Vq — ftm accompanied by be- 
ginning to move the left barrier (optical piston) with a 
velocity v p = 0.01 to the right and phase imprinting the 
matching velocity v p onto the condensate. As our inves- 
tigations here require large reservoir extensions and long 
simulation times, we resort to a one dimensional model 
in this and the following subsection, employing the po- 
tential V(f = 0,x,t) in Eq. (|2"TT) . 

Fig. [3 highlights the differences between these meth- 
ods. Consider panel (a). When the barrier at i r — 80 is 
suddenly reduced, the interaction energy Uh exceeds Vq 
and the BEC begins to stream out. A wave within the 
reservoir, resulting from the sudden change of the poten- 
tial, propagates left, hits the barrier there and reflects. 
In contrast the initiation with the piston (panel b) shows 
less perturbation of the reservoir. Panel (c) shows a Mach 
number profile near the hump potential which represents 
either case well. Consider the analogue Hawking tem- 
perature as a function of time shown in panel (d). For 
this and similar graphs, we numerically determine the 
instantaneous location of the horizon and evaluate the 
dimensionless horizon temperature Eq. (|12p there. 

During the timespan covered in panel (d), the speed 
of sound at the horizon decreases by 40% if no piston is 
used. Nonetheless, as predicted by Eqs. and (TIT))) , the 
analogue Hawking temperature remains fairly constant, 
as the decrease in speed of sound is mostly balanced by 
a corresponding increase in the Mach number gradient. 
However, step like decreases in temperature occur. The 
steps coincide with successive arrivals of the flow pertur- 
bation, visible in panel (a), which is repeatedly reflected 
from the left and right barriers. One has to keep in mind 
that Eqs. ([5]) and (TIT))) were derived neglecting the quan- 
tum pressure term Q(x) and assuming strictly a station- 
ary state. 
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FIG. 2: Initial time evolution of the condensate in the reser- 
voir after initiation of the outflow. No three-body loss, (a) 
Density (solid) and V(x)/Uid (dashed) for outflow with no 
piston. p, 1D = 0.61, i bg = 0.64, & = 3.8, a = 14 ~ 226 s , 
Vo = 0.4, x r = 80, xi — —110, d = 0.55-. The time samples 
are t~0, 60, 210, and are shifted up for clarity. The zero of 
the common axes for h and V for later time samples can be 
deduced from V"(0) =0. f * = 0.008. (b) The same as in 
(a), but with an optical piston. Here Vq = 0.61, £(, 9 = 0.64, 
£h = 4.66, d = cr, all other parameters are as in panel (a). 
T* = 0.005. (c) Sample Mach number profile (solid, thick) 
for the case with piston. The thin black line indicates M = 1. 
(d) Comparison of analogue Hawking temperatures with (o) 
and without (x) piston. 



It is evident in Fig. [2] that the density at the horizon 
is significantly less than the bulk value. The value of the 
healing length at the horizon £/, indicates that the hori- 
zon in fact slightly violates our demand Eq. (fT3| . This 
does not affect the conclusions of the present section. 
The big £h has two important consequences: Firstly, the 
quantum pressure term Q{x) begins to have an effect 
on the condensate background flow. The effect is small: 
for Fig. O Q(xh) initially makes up only 2% of the en- 
ergy components in Bernoulli's equation Secondly, 
the non-phononic part of Bogoliubov's dispersion rela- 
tion has an effect on phonon propagation, quantified in 
Ref. [l5j]. This might have to be tolerated for an outflow 
experiment for the following reason: The bulk conden- 
sate in Fig. [5] is as dense as allowed by loss. Thus the 
only way to increase the density at the horizon, and hence 
the temperature, is by reducing Vq in Eq. (|24p . However, 
if Vq is reduced, the bulk perturbation created during 
initiation becomes much larger. 
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FIG. 3: (Color online) Time evolution of analogue Hawking 
temperature Th for Jx\d = 0.61, = 0.64, £ h — 2.8, a — 
14 ~ 22^j, s , Vb = 0.5. Initially for these scenarios T* = 
0.014. (a) Left end barrier acts as piston with v p = 0.01. 
(x-black) K 3 , 1D = 0. (o-red) K 3 , 1D = 1.4 x 10~ 7 . (b) Left 
end barrier does not move, (x-black) Kz,id = 0. (o-red) 
Kz.id = 1-4 X 10~ 7 . t = 5000 corresponds to t = 117 ms. The 
three-body-loss constants are for helium with c^x = 68 00 x 27r 
Hz transverse confinement (table [I]). 



D. Effect of Three-Body-Loss 

In this section, we show that if inelastic atom colli- 
sions are taken into account, the piston is no longer nec- 
essary. Without the piston, and without three-body loss, 
the reservoir density hb g decreases in time, because the 
reservoir is emptied by the outflow. However, for a large 
reservoir the fractional atom loss due to the outflow, and 
the resulting density decrease, will be small. Whether or 
not a piston is used, we expect strong three-body losses 
in the high density regime, which we are forced to enter 
to reach high Hawking temperatures. These losses will 
dominate density reduction. 

In this section, we employ parameters for helium, as- 
sume uj± = 6800 x 2ir Hz and a bulk density correspond- 
ing to n = 3.6 x 10 19 m -3 , roughly as in table [IT] We 
use a large reservoir with xi = —2270, x r — 80, d = a/2. 
Further parameters are given in the caption of Fig. [3] 

We find a significant reduction in the analogue Hawk- 
ing temperature on the time scale of 50 ms. While these 
simulations begin from a state safely in the hydrody- 
namic regime, within the limits derived in section IIII1 
it does not remain so during the time-evolution: Three- 
body loss reduces the bulk density and hence also the 
density at the horizon. As a consequence the local heal- 
ing length near the horizon approaches the width of the 
hump potential, hence the flow ceases to obey the hydro- 
dynamic Eq. (0. 

We note that the time evolution of the analogue Hawk- 
ing temperature shows no advantage of the piston setup 
over free outflow, see Fig. [3] The bulk interaction energy 
of the initial state condensate Uiofibg = [Aid = 0.61 is 
significantly larger than the exit barrier height Vb = 0.5, 
giving some room for a reduction of condensate density 
due to three-body losses without an interruption of out- 
flow. 



E. Three-dimensional Horizon Simulations 

Finally, we report on a 3D simulation including loss, 
which shows the feasibility of the proposed experiment. 
The parameters were chosen to represent the case with 
the highest analogue Hawking temperature of table HT1 a 
sodium condensate with uj± = 6800 x 2-7T Hz (a osc = 
2.5 x 10~ 7 m), resulting in h^rj = 30, which corre- 
sponds to = 1.8 x 10 21 m~ 3 . This gives jl^D = 4.3, 
which places the simulation between the Q1D and TTF 
cases. We choose a scenario that is not in the quasi one- 
dimensional regime, to show that even if transverse dy- 
namics are possible, a stable horizon is obtained. Further 
parameters were Vb = 2.8, a — 6.8. The healing lengths 
are £ bg = 0.36 and £ h = 1.4. T* x 50 = 2.8. Our simu- 
lated BEC cloud is about 4 units (1 /an) in diameter and 
70 units (17.5 jum) long. Following the conclusion of the 
previous section, we do not use an optical piston. 

We find stable outflow with Tjj = 0.07 (evaluated 
at r = 0), in good agreement with the expectation of 
Th = 0.068 from Eq. (f2"2"| . For the assumed transverse 
confinement this equals 23 nK, somewhat exceeding the 
estimate of the achievable temperature given in table HT1 
It decays slightly over the timespan of t = 600 (14 ms) 
of the simulation, largely due to the condensate outflow. 
Computational limitations required a small spatial do- 
main (17.5 /xm), containing only 30000 atoms. In an 
experiment the length of the reservoir could be much 
larger, reducing the temperature reduction due to con- 
densate outflow as shown in section HV Dl 

A comparison of the atom number evolution with and 
without three-body loss separates the three-body loss 
contribution to the reduction of the atom number from 
that due to outflow. We find that 5% of our initial atoms 
are lost due to three-body recombination. 

Our simulation suggests that the Hawking tempera- 
tures in tablellllcan be realized under realistic conditions, 
even in a not entirely one-dimensional situation. Fig. 0] 
shows the two dimensional structure of the BEC in the 
rx plane as well as the radial variation of density, Mach 
number and the approximate horizon temperature [33j 
T p = \M'\c/(2tt) for a slice of constant x — Xh- This 
slice does not coincide with the sonic horizon, whose x 
co-ordinate is a function of r. Finding the exact horizon 
location is difficult. The violet line in Fig.0](a) is the in- 
tersection of the rx plane with the surface where M = 1, 
which is the boundary of the ergoregion [34| . On this sur- 
face, the flow velocity equals the speed of sound, but since 
the flow is not perpendicular to the surface, the velocity's 
normal component does not. The true horizon slices the 
rx plane somewhere between the violet and white lines. 
Due to the spatial proximity of the slice shown in Fig. Q] 
and the true horizon, we expect the displayed behaviour 
of T p to give a very good indication of the transverse 
temperature structure at the true horizon. At r = the 
graph shows the true horizon temperature Th- 

The Mach number increases away from the trap axis, 
r = 0. This is largely because the density decreases while 
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the velocity is approximately independent of r. T p stays 
roughly constant on the shown slice. This qualitative 
behavior is found also for simulations with different pa- 
rameters. 

For the shown scenario, the condensate healing length 
is £ = 0.36, which is less than the transverse extension of 
the BEC, about R ~ 2. The transverse structure could 
thus be resolved by an excitation with wavelength A ~ 
l/£, which might affect the excitation's behavior. In a 
Q1D scenario, the horizon also has a curved shape, but 
in that case £ ^> R. Hence the transverse dimension is 
frozen out for phonons and the radial structure of the 
horizon should be irrelevant. 

The horizon is stable in three dimensions, whether or 
not the confinement is tight, and for a condensate beyond 
the Q1D regime its transverse shape is non-trivial. 
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FIG. 4: (Color online) Transverse structure of the BEC near 
the horizon. The full length of the condensate is about 70 
units. The hump potential is centered at x r = —4 with 
width a = 6.8 [35|] . The snapshot is taken at t — 525 (12 
ms after initation). Other parameters: xi = — 76, d = 0.7<r. 
(a) Colormap of density in rx plane. The density is high- 
est on the left (shaded red) and monotonically decreases to 
the right (shaded blue). Superimposed are contours of equal 
Mach number (red) at M=0.5, 2, 4. We marked M = 1 in 
violet, (b) Radial structure of Mach number M (dot-dashed), 
density n (dashed) and proto analogue Hawking temperature 
T p (solid, see text for definition) at x = Xh — —3.54. This 
slice through the condensate is indicated by the white straight 
line in (a). 
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